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Abstract: 

Introduction: In a real-world scenario, problems are identified, assumptions 
are made, data are collected and models are proposed. The underlying math- 
ematical structure of the models is analyzed to appraise the sensitivity of the 
conclusions when the assumptions are not precisely met. 

Aims: In this work, we present the formulation of a modified Lotka-Volterra 
Predator-Prey Model to include terms which represent the degrees of inter- 
nal competition of both predators and prey for their limited resources and a 
simulation of the stability behaviour of the resulting system. 

Materials and Methods: A computer simulation of the behaviour - stability 
of the solution of the model - of the system under different conditions is also 
presented. The critical points, eigenvalues and phase portrait of the system 
in terms of the model parameters are also obtained. 

Results: The model interpretation and applicability are also discussed. 
Conclusion: The trajectories of the proposed model are not periodic and 
tend to equilibrium level with time, thus, making the system to be asymptoti- 
cally stable. This behaviour is in line with the observable natural predator-prey 
system. 

Keywords: Lotka-Volterra, Predator-Prey Model, Simulation, Critical Points, 
Eigenvalues, Stability 


All co-authors agreed to have their names listed as authors. 


This is an open access article under the terms of the Creative Commons Attribution License, which permits use, 
distribution and reproduction in any medium, provided the original work is properly cited. 
© 2018 The Authors. Journal of Research and Reviews in Science — JRRS, A Publication of Lagos State University 


JRRS 


https:sciencejournal.lasucomputerscience.com 


LASU Journal of Research and Review in Science 


1 INTRODUCTION 


Interest in how populations tend to grow gained attention in 
the 18th century when the work of the author in [f] was pub- 
lished. The author in his work [i], put forth an exponential 
growth model for human population and concluded that even- 
tually the population would exceed the capacity to grow an 
adequate food supply. 


However, the assumptions of the model left out factors (envi- 
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2 MATERIAL AND METHODS 


2.1 Identification of Problem 


The author in [7] noted that as more individuals are produced 
than can possibly survive, there must in every case be a 
struggle for existence. This is true as the individuals are faced 
with the need to compete for their limited resources either 
with species of same kind or distinct species. In this work, 


ronment and interaction between population of different species We take a closer look at the situation between the predator 


that are important to population growth and hence it could not 
stand the test of time as it was proven to be inaccurate for 
technologically developed countries. 


Interestingly, the work [1] stimulated interest in the study of 
the ecological problem of a predator population interacting 
with the prey one which has been of great interest in many 
areas of science. Early attempt to study this problem us- 
ing mathematical equations is attributed to the independent 
works of [2] and [8] 4]. The authors independently produced 
the equations that gave the model of this problem and dis- 
covered that, under simple hypotheses, periodic fluctuations 
of the populations occur. Hence, the mathematical equations 
which model the struggle for existence of two species (preda- 
tor and prey) bear the names of the pioneer scientists: Lotka 
and Volterra. The factors that were left out in the model pro- 
posed by the author in [1] formed the basis for the formulation 
of the Lotka-Volterra model. The assumptions about the en- 
vironment and evolution of the predator and prey population 
for the Lotka-Volterra model are: 


e The prey population have an unlimited food supply at 
all times. 


The prey population x(t) increases at arate dx = ax(t)dt 
(proportional to the number of prey) but is simultane- 
ously destroyed by predators at arate dx = —ba(t)y(t)dt 
(proportional to the product of the numbers of predator 
and prey). 


The predator population y(t) decreases at a rate dy = 
—cy(t)dt (proportional to the number of predators), but 
increases at a rate dy = dx(t)y(t)dt (again proportional 
to the product of the numbers of predator and prey). 


The first and simplest predator-prey model ({2], [8], [5], [6]) is 
given by the nonlinear system of first-order ordinary differen- 
tial equations 


= = ax(t) — bæ(t)y(t) (1) 
ou = —cy(t) + dx(t)y(t), 


In this work, we present a simulation of the stability behaviour 
of a modified version of the predator-prey model. 
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and prey. The predators eat the prey and the prey live on the 
abundantly available natural food in the environment. If the 
predators eat so many of the prey that the prey cease to be 
abundant, the food supply of the predators is greatly reduced, 
then the predator will starve. A dwindling in the population of 
the predator will result in a significant increase in the prey’s 
population. As the prey population increases, the food sup- 
ply of the predator grows and so does the predators’ popula- 
tion. Also, more predators are eating increasingly more prey 
again. In this cyclic situation, one asks whether the cycle 
continues indefinitely or does one of the species eventually 
go into extinction. The ability to answer the above question 
finds relevance in many areas of science. 


2.2 Model Formulation and Assumptions 


It was observed that internal competitions between predator 
and prey for their limited resources were not considered as 
major factors in the models earlier proposed. In this work, we 
considered a modified Lotka-Volterra equation by introducing 
the degree of internal competition of the predator and prey 
for their limited resources. We assume that the degree of 
internal competition of the predator and prey are proportional 
to the square of the predator and prey populations, i.e. fy?(t) 
and ex?(t) respectively. Hence, the modified Lotka-Volterra 
predator-prey model presented in this work is given as 


< = ax(t) — ba(t)y(t) — €0(t) > 
= —cy(t) + de(t)y(t) — fy?(t), 


where a,b,c,d, e, f are positive constants. 


2.3 Stability Property of Model Solution 


To analyse the behaviour of the system under consideration, 
some properties of the system for varying parameters of the 
model need to be investigated. Here, we present the critical 
points, eigenvalues, phase portrait and stability of the model 
under consideration. 


2.3.1 Critical Points of Model 


To obtain the critical points of the model (2), we set dx/dt = 
f(x,y) = 0 and dy/dt = g(x,y) = 0 and solve. Unlike a lin- 
ear system, a nonlinear system could have none, one, two, 
three, or any number of critical points. Since there might be 
multiple critical points present on the phase portrait, each tra- 
jectory could be influenced by more than one critical point. 
This results in a much more chaotic appearance of the phase 
portrait. Consequently, the type and stability of each critical 
point need to be determined locally on a case-by- case basis. 
Considering the system 


dx 
dt 


oY = oyli) + detu) — fy?) = 0, 

Since the solutions x(t) and y(t) are not easily obtained by 
mere inspection, we proceed with the use of Mathematica”’s 
built-in command for solving equations. Solving the system 
above for x(t) and y(t) results in four critical points of the 
system. Hence, the critical point of the model @) are 


em (0-9) 


—af — be 
(an) = ( bd+ef’ 


(x3, y3) = (0,0) 
Cnu) = ($0) 


Interestingly, three of the critical points of the system are de- 
pendent on the values of the parameters a,b,c,d,e,f. 


= ax(t) — bz(t)y(t) — ex?(t) = 0 
(3) 





ce — ad 
a (4) 


2.3.2 Eigenvalues of Model 


A 


is called the Jacobian matrix of the system. It contains the 
first order partial derivatives of f and g evaluated at the point 
a — by — 2ex 


(x,y). Now 
ia ( ) © 


Evaluating 6) respectively at the four critical points @) results 


In 
a+ be 
h(s) = ( af 2: (7) 


The matrix 


On f(x,y) 


PEE BN: Goa a) 


Oy f(x,y) ) (5) 


—br 
—c + dz — 2 fy 








e(af+bc) b(af+bc) 
bd bd 
J2 (x, y) = ( Pica cof ad (8) 
bd-tef bdFef 
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nen = (4 (9) 


ab 
nen = (F ate) (10) 


The eigenvalues of (7), (8), @) and (H0) respectively are 




















AL = cC 
_ af +bc (11) 
My a 
rA2 = T (—bce — adf — aef + cef— 
y ((adf + aef + bce — cef)? — 
4 (abd? f + a°def? + abcd” — ace? f? — bc*de — be?e? f))) 
u2 = T (—bce — adf — aef + cef + 
y ((adf + aef + bce — cef)? — 
4 (a°bd? f + a?def? + abcd? — ace? f? — bde bee? f))) 
(12) 
A3 a (13) 
3 = =c 
Ag = —a 
ad — ce (14) 





2.4 Stability of Model 


We can use the eigenvalues of the Jacobian matrix to decide 
the type and stability of the critical point (x,y). Let A and u be 
the two eigenvalues. Assume that à and yu are different from 
0. 


2.4.1 Phase Portrait and Field Directions of the Model 


The best way to describe the vector field of the system is to 
simulate the plot with varying values of the parameter. Since 
a static plot of the vector field for the model can only be ob- 
tained by fixing the values of the parameter, we present some 
of the possible scenarios with three case: 
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Table 1: Conditions on eigenvalues to ensure different stability 





Eigenvalues Linear System 











Nonlinear System 
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Figure 1: Phase portrait and population plots for @) with a = 2,b=1.1,c=1,d=0.9,e= f =0.1 


3 RESULTS 


3.1 Model Interpretation and Conclusion 


Figures [I], 2 and | show that the stability at the various criti- 
cal depends on the values of the model’s parameters. Based 
on the introduction of the degree of internal competition be- 
tween the predator and prey, we see from Figures [1], 2 and B}, 
that the trajectories of the model are not periodic and tend to 
equilibrium level with time. Hence, the system B) is asymp- 
totically stable. This behaviour is in line with the observable 
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natural predator-prey system. 


4 CONCLUSION 


We have been able to modify the classical Lotka—Volterra 
equation by introducing terms which represent degrees of in- 
ternal competition between the predator and prey. Unlike the 
classical model, this modification has been shown to make 
the system asymptotically stable. This behaviour conforms 
to the observable natural predator-prey system. For future 
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Figure 2: Phase portrait and population plots for B) with a = 0.9,b = 0.1,c = 0.1,d = 0.1,e = f = 0.1 
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Figure 3: Phase portrait and population plots for @) with a = 0.9,b = 0.1,c = 0.1,d =0.1,e = 2.5, f = 1.9 


work, more variables can be included in order to get a better AUTHORS’ CONTRIBUTIONS 
representation of the system. 
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